Correlations of Rydberg excitations in an ultra-cold gas after an echo sequence 
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We show that Rydberg states in an ultra-cold gas can be excited with strongly preferred nearest- 
neighbor distance if densities are well below saturation. The scheme makes use of an echo sequence 
in which the first half of a laser pulse excites Rydberg states while the second half returns atoms to 
the ground state, as in the experiment of Raitzsch et al. [Phys. Rev. Lett. 100 (2008) 013002]. Near 
to the end of the echo sequence, almost any remaining Rydberg atom is separated from its next- 
neighbor Rydberg atom by a distance slightly larger than the instantaneous blockade radius half-way 
through the pulse. These correlations lead to large deviations of the atom counting statistics from 
a Poissonian distribution. Our results are based on the exact quantum evolution of samples with 
small numbers of atoms. We finally demonstrate the utility of the tj-expansion for the approximate 
description of correlation dynamics through an echo sequence. 

PACS numbers: 32.80.Ee, 32.80.Rm, 34.20.Cf 



I. INTRODUCTION 

When atoms within an ultra-cold gas are excited to 
Rydberg levels, they experience long-range interactions 
which can block further excitations, leading to a strongly 
correlated many-body state. This effect might be useful 
for quantum information H as well as for fundamental 
studies of maw-body physics and was observed in several 
experiments 0, i, H IE IE 0, H i, EE El • The coherence 
of the excitation process in a bulk gas has been exper- 
imentally demonstrated using an echo-technique [1, 0: 
After exciting atoms to Rydberg states, it was possible 
to de-excite them following a 7r-phase shift of the exci- 
tation laser. The basic scheme is illustrated in Fig. [1] 
After the echo sequence, a certain fraction of the atoms 
remains in the excited state, owing to effects of the in- 
teraction. This has been modeled theoretically using the 
super-atom approach [12|. 

Neither experiment nor theory has considered the dy- 
namics of Rydberg-Rydberg correlations during such an 
echo-sequence. Here, we show that strong correlations of 
atoms separated by a characteristic distance rg are in- 
duced in its course. This distance ro is slightly larger 
than the instantaneous blockade radius at the moment 
when the laser phase is flipped. At the instantaneous 
blockade radius rb(t) the Rydberg state density-density 
correlation function drops sharply to zero. Initially rb(t) 
equals zero, later growing towards its saturation value 
7-&0 ~ (Ce/H) 1 / 6 , as the longer duration of the pulse al- 
lows an increasingly finer energy resolution. 

After most Rydberg atoms that were excited in the 
first half of the pulse returned to the ground state, the 
majority of the remainder has their nearest excited neigh- 



bor at a distance between ro and 1.5ro. The strength of 
this correlation signal is proportional to p -2 , where p is 
the atomic density. 




FIG. 1: Illustration of an echo excitation sequence and evo- 
lution of the Rydberg fraction. For a sequence of duration 
T, we use the fractional time r = t/T. (a) Rabi frequency 
Qw(t) [see section III A] during the sequence, (b) Fraction 
of excited state atoms, /, for a density poVbo = 26.2 (solid) 
from a numerical solution of Schrodinger's equation. Here, 
Vi,o = 47rr^ /3 is a simple estimate of the saturated block- 
ade volume. We also show the non-interacting result for the 
excited fraction (dashed). The symbols (•) mark times for 
which we display correlations in Fig. [2] 
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Our results imply that the echo technique can be used 
to manipulate the nearest-neighbor distribution in the 
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Rydberg fraction of the gas. Such manipulations could 
be used for example to initiate dynamics due to dipole- 
dipole forces [l3[ from a well specified non-equilibrium 
state. 

Similar correlations between Rydberg atoms are cre- 
ated if a two-step excitation scheme with a strongly de- 
caying intermediate state is used flil . [Til ] . In that case 
they are due to the Autler-Townes splitting of the inter- 
mediate level. 

The peculiar nature of the nearest neighbor distribu- 
tion function of Rydberg atoms after an echo pulse could 
be seen experimentally in the kinetic energy spectrum 
of ions after field ionization of the remaining Rydberg 
fraction. Alternatively one could measure deviations of 
the atom counting statistics from a Poissonian distribu- 
tion [l|| . Both these suggestions are discussed further in 
section IIVI 

Our results are based on solutions of the many-particle 
Schrodinger equation, which we also use to benchmark 
the recently proposed w-expansion [l?], QjJ ■ Both meth- 
ods are briefly described in section [TTJ The ensuing cor- 
relation dynamics are presented in section IIII1 Possible 
ways to detect the pairing effect are discussed in sec- 
tion [IV] In section [V] we take a closer look at the effects 
of density variations on our results and in section IVII we 
directly compare many-body quantum simulations with 
the w-expansion. Finally, we conclude in section fVIIl 



takes the form: 



i i,j;i>j 

- (nw{t)c n ,n + fi*w*(i)c n , ti ) , (2) 



where c n .j; (c n ,fi) is the coefficient of the state that is 
reached from |n) by lowering (raising) the i'th atom if 
this is possible, and otherwise. 

For the numerical solution we convert Eq. (|2|) to a ro- 
tating frame for the variables 



(3) 



To render simulations of ensembles with hundreds of 
atoms possible, we remove all states whose interac- 
tion energy is larger than some cut-off E cu t from the 
Hilbertspace (see eg. [l9, 20]). We further limit the num- 
ber of simultaneously excited atoms to ^n,- < -/Vmax- 
All numerical results are checked for convergence with 
respect to variations of N max and E cut . 



II. METHODS 

A. Exact quantum dynamics 

We consider a system of Nq atoms at fixed locations xi , 
described by the following Hamiltonian in atomic units: 

j No 

i=l 

N No 
+ A E< + £ K H<e4e- (!) 

Each atom can either be in its ground state \g) or an 
excited state |e), the Rydberg state. The atoms in |e) 
experience long range interactions Ky = — C§j\xi — xj\ e . 
Laser induced transitions between the levels occur with 
Rabi-frequency O G R, detuning A and the temporal 
profile of the pulse w(t). The operators a\ h act as o % ab — 
\a)(b\ on the subspace of atom i and as unity on others. 
We use atomic units unless otherwise indicated. 

We expand the many-body wave function as = 
S n Cn l n )' wnere n is a vector with N elements rij S 
{0,1} which describe whether atom j is in the ground 
(0) or excited (1) state. Schrodinger's equation (SE) then 



B. The ^-expansion 

In the experiment ||, the sign of the Rabi-coupling 
was flipped well before the excitation of Rydberg atoms 
reached the saturation limit imposed by the blockade ef- 
fect. In such a situation one could attempt to solve the 
quantum dynamics of the gas by a series expansion in 
lu = fit [TiT ]. where to — 0.25 for the experimental condi- 
tions. 

From Eq. {!]), we can derive equations of motion for 
the operators er^ b . We expand these operators in a se- 
ries a l ab = J2n ujn,J ab^ ■ ft i s possible to obtain recursive 
equations expressing by er^,™' ) for m < n. In our 

case the lowest order approximations of the a l ab already 
offer interesting insight into correlation dynamics. 

In the following we are interested in an echo- type pulse. 
Assuming a pulse-length T, we rescale our variables 

r = t/T, uj = QT, % = Kij T, /(r) = w{t/T). (4) 

The echo pulse shown in Fig. [T]is then given by /(t) = 
6(1/2 - r) - 9(r - 1/2), where 9(f) is the Heaviside 
function. We further define F(t) = J Q dr'/(r') = 
t9(1/2 - t) + [1 - t]9(t - 1/2). 

According to the w-expansion the leading- (LO) and 
next-to-leading- (NLO) order expressions for the excited 
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case 


(i) (ii) 


(iii) 


(iv) 


N 


20 84 


125 


50 


L/r b0 


2.5 2.38 


2.38 


2.38 


Pg{r) 


po po 


Po 


po exp [-2r 2 /a 2 ] 


PoVbo 


5.4 26.2 


39 


26.2 


po/10 10 [cm- 3 ] 


1.3 6.3 


9.4 


6.3 



TABLE I: Parameters for the scenarios modeled in this article. 
In cases (i)-(iii) the atoms are homogeneously distributed in a 
cubic box of side-length L with periodic boundary conditions 
[2lj . Case (iv) is a gaussian cloud as shown, with a — 1.6rf,o- 
This ensures the same peak-density as the homogenous den- 
sity of case (ii). The relevant dimensionless measure of the 
density is poHo- For the parameters described in the text this 
would correspond to the values in the last row. 



compare their predictions with direct numerical simula- 
tions in section IVI1 A more detailed description of the 
w-expansion can be found in (l7j . 

To estimate which order of the expansion in lo is re- 
quired for a given scenario, we can consider the excited 
fraction of a fully blockaded sample of Nf, atoms assum- 
ing a square pulse, 

P e {t) = ^sin 2 [v/AW2] =^pW(fit)", (12) 

n 

which is correctly reproduced by the expansion in the 
limit of infinite interactions fl7| . Roughly knowing the 
expected Nb for a system and the maximal fii, we can 
estimate how many terms of series expansion Eq. fj 1 2[) 
are required. This will be used in section IVll 



state fraction P e (r) are given by: 
P^°\r)=^{a^)+^{a^) 

hi = — ^ Re 1^—^— dnfin) F(n) 

742 = lE R 4 [ T dnf(n)[F(r) - 2F( n )} 



(5) 

(6) 
(7) 



< J dT 2 r(T 2 )F*(T 2 )(e* 



(Tl-T 2 )/8i. 



1 



(8) 



The Rydberg-Rydberg correlation function is defined by 



{oU(r)){aie{r)) 
with leading order approximation in the w-expansion 

4\^dn e^7(Ti)F(Ti) " 



(9) 



\F(r) 



(10) 



Calculating higher order corrections to these quanti- 
ties, although possible in principle, is more tedious than 
justified. The expressions above can be explicitly evalu- 
ated for a homogeneous system with van-der-Waals inter- 
action Kij — — Ce 1 1 xi — Xj ; | 6 . The sum in I± 2 is replaced by 
Y^ijtj ~ * f d 3 xp(x), corresponding to an ensemble aver- 
age. Throughout the paper we write p for the full density 
profile and po for the peak density. If we evaluate Eq. ([6]) 
at t = T for a homogeous system, we can obtain the echo 
signal, i.e. the final Rydberg fraction 



f e = P e (T) = 



2835 



(8^2- 9) p . (11) 



Other analytical expressions that can be obtained are 
not illuminating and hence omitted here. However, we 



III. CORRELATION DYNAMICS 

For A^o atoms, homogenously distributed in a cube of 
volume L 3 , we can see that the physics of our problem 
is governed by two parameters, namely lu — CIT and 
U = C 6 T/L 6 . Throughout this paper we usually em- 
ploy u> = 0.25 and U = —0.0014. This corresponds, for 
example, to a 87 Rb gas in which states with principal 
quantum number riRyd = 41 are excited via a transition 
with Rabi- frequency £1 = 2tt x 0.1 MHz during a time of 
400 ns. The side- length of the box would be L — 11pm. 
For these parameters we then vary the number of atoms 
as listed in table U leading to densities of the order of 
5 x 10 10 cm~ 3 . Compared to the experiment [8[, these 
parameters lead to a substantially weaker blockade with 
maximally about Nf, = 40 atoms per blockade-sphere. 
This facilitates both, the physics of interest here and the 
numerical simulation. We will use the saturated two- 
atom blockade radius r^o — (Cq/Q) 1 ^ 6 as length scale 
and express densities using the simple estimate for the 
blockade volume Vj, = 47rr 3 /3. 

In an echo-sequence, we first excite Rydberg states 
with Rabi- frequency fl for some duration r/2, followed 
by a 7r-phase shift of the laser and hence Rabi coupling 
— for further time r/2. These parameters are defined 
in section Til Bl For the full quantum dynamics as in sec- 
t ion III Al we randomly distribute No atoms, solve the SE 
and obtain one correlation function defined by Eq. §9§ 
for each pair of atoms. These are binned according to 
the separation r of those atoms. To minimize finite size 
effects, we only consider atoms (pairs) inside a central 
cube of side length L /2 within our simulation volume to 
calculate N e (g [2) ). 

In order to obtain a spatial correlation function, we 
first calculate g^ 2 '(i,j) according to Eq. ((9]) for each of 
these atom pairs. We then determine 



(13) 



where Y^i j denotes a double sum over all atoms that 
fulfill \xi — Xj\ G [r, r + Ar) for a bin size Ar, and N(r) 
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tance r. Let us write the quantum state of this pair as 



FIG. 2: (color online) Spatial atomic correlations during an 
echo sequence with poVbo = 5.4. (a) Pair correlation function 
g^ 2 \r) during the excitation phase at r = 0.05 (black), r = 
0.37 (blue) and r = 0.49 (red). Below the cutoff length r cu t, 
is set to zero [21JJ . (b) The same during the de-excitation 
phase at t = 0.6 (black), r = 0.7 (blue) and r = 0.75 (red), 
(c) Final shape of the pair correlation function after the pulse 
(r = 1). (d) Spatial maximum of pair correlations. The 
location of the time-samples in (a-c) relative to the pulse can 
be seen clearly in Fig. [T] 



is the number of atom pairs that fall into each respective 
bin. 

Thus averaging over the whole sample and further over 
a large number of realizations of the spatial atomic dis- 
tribution [221 ] . we obtain a spatial correlation function 
g^(r). Our results from numerical solutions [21| of 
Eq. (2j) are shown in Fig. Q] and Fig. [2] 

Soon after initiating a transfer of atoms to the Ryd- 
berg state, the probability to find a pair of excited atoms 
closer than the sharply defined instantaneous blockade 
radius r&(t) is essentially zero. This is due to strong van- 
der-Waals interactions of those pairs, shifting doubly ex- 
cited states out of resonance. The radius rb(t) grows in 
time, since the longer pulse duration allows an increas- 
ingly finer energy resolution. The Rydberg atom number 
in the echo-sequence does not reach the saturation limit, 
hence Tf,(t) does not reach its equilibrium extension (rbo) 
but continues to grow throughout the excitation period 
[Fig. [2j panel (a)]. Following the flip of the phase of the 
Rabi-coupling at r = 0.5, </ 2 ' develops a dominant peak 
just outside the blockade radius [panel (&)]. The height 
of this peak increases as more and more atoms are trans- 
ferred back to the ground state [panel (d)]. Just before 
the end of the pulse, for densities as low as in Fig. a 
dip appears in the correlation peak like that shown in 
panel (c). This feature starts to occur at the time when 
the probability for a single excitation in the system drops 
below that for double excitation. 

The underlying physics requires only binary interac- 
tion: Consider a pair of atoms separated by some dis- 



I*) = c gg\9i92) + c eg \eig 2 } + c ge \gie 2 } + c ee |eie 2 ) 



(14) 



Only in the state \e1e2} where both atoms are excited to 
a Rydberg level do they experience any interaction. The 
correlation function in state (Tl4"l) is 



5 (2) (1,2) 



<-y 



'■)(\c 



.</<' 



+ C e 



(15) 



Now consider g( 2 > after the pulse for three different 
atomic separations: (i) If the atoms are very close 
r < 7"b(£) double excitation can be considered fully sup- 
pressed. Hence |c ee | = and </ 2 ^ = 0. (ii) If the atoms 
are very far apart the interaction can have no effect. We 
know then — 1. From Eq. (fTS"]) this can be under- 
stood since without interactions and for small excited 



fractions \c e 



-egl 



and 



-eg I 



(iii) 



In the intermediate range we can neither neglect dou- 
ble excitations nor interaction. Consider the very end of 
the pulse: The amplitudes c eg and c ge have returned to 
their initial value of zero after the echo-pulse. In contrast 
|c ee | is non-zero due to the dephasing. We then see that 
the atomic correlation function scales as g^ = l/|c ee | 2 , 
which is larger than one. 

As shown in Fig. [2j the distances r with pairing cor- 
relations during the second half of the pulse are those 
where correlations change, as r increases, from blockaded 
(g^ = 0) to uncorrelated (g^ = 1) during the first half 
of the pulse. We will call the spherical shell around each 
atom where neighbors have these distances the partial 
blockade shell. 

We find that the position of the maximum of the cor- 
relation function during the second half of the pulse de- 
pends only weakly on time. 



IV. SIGNATURES OF PAIRING 

In this section we quantify to what extent the corre- 
lation dynamics presented in the previous section allows 
control over the nearest neighbor distribution in a Ryd- 
berg gas. We further discuss observables that are easier 
accessible experimentally than the density-density corre- 
lations themselves. 

We consider the fraction R of excited atoms separated 
from their nearest neighbor by a distance in the interval 
[ r o. r o + d], chosen to contain the peaks in Fig. [2] (a). This 
is given by 



R 



En # o|Cn| 2 ;(n) 
En^O l C n| 2 



(16) 
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FIG. 3: Visibility of the pairing effect for a density with 
pVbo = 26.2. (a) Mandel-Q parameter, (b) Fraction R of rem- 
nant atoms within the preferred distance peak, using Eq. (11611 . 
(c) Total number of paired, excited atoms, extrapolated [23| 
as if there was a total number of N g = 1 x 10 7 atoms initially. 



(ii) The Q-factor captures only integrated properties of 
the correlation function. To obtain information about the 
spatial shape of atomic correlations after an echo pulse, 
the paired Rydberg atoms could be field ionized. The 
potential energy of the ions due to Coulomb repulsion 
Seoul ~ e 2 /ro47reo would subsequently be converted into 
kinetic energy. Since only a small number of atoms oc- 
cupies Rydberg states, one could hope that the Rydberg 
fraction is sufficiently dilute for Coulomb-scattering to 
leave the initial kinetic energy distribution essentially un- 
changed. The peak in the nearest-neighbor distribution 
function ~ r 2 g^ 2 \r) then translates into a well visible 
maximum of the measured kinetic energy spectrum of 
the ions. For pVbo = 26.2 the energy at the maximum is 
0.5 meV. 



where /(n) is the fraction of excited atoms in the many- 
body basis state |n) with at least one excited neighbor 
in the chosen interval. 

For a situation as in Fig. [2] we choose r = 0.5r{, 
and d — 0.25r(,o- We show in Fig. [3] (b) that almost 
all atoms are paired up before the end of the pulse. In 
the initial phase, R remains nearly zero owing to the 
predominance of single excited atoms without neighbor. 
In panel c we show the total number of paired, excited 
atoms, obtained by multiplying the paired fraction R(t) 
with the Rydberg- number N e (t). We then rescale the 
obtained number corresponding to a total initial num- 
ber of 10 7 atoms for illustrative purposes. It may seem 
counter-intuitive that the number of paired atoms rises 
even during the second half of the pulse when atoms are 
predominantly deexcited. Note however, that whether or 
not the Rabi-coupling causes excitations or de-excitations 
depends on the relative populations and phases. The 
continuous increase of RN e (t) is again due to interaction 
induced decoherence. 

Finally we describe two possibilities to experimentally 
detect the correlation dynamics: (i) Via the effect of cor- 
relations on the number uncertainty in the Rydberg frac- 
tion [l6| and (ii) via field ionization of the paired Rydberg 
atoms. 

(i) Strong deviations from an uncorrelated state have 
been shown to affect the number statistics of the excited 



state fraction 
Q parameter [24 1 




This is well captured in the Mandel- 



d 3 r Pe (r)[gW(r)-l], (17) 



which can be experimentally determined from the Ryd- 
berg atom counting statistics. We show the Q-factor in 
Fig. [3] (a). In comparison with the evolution of R, we 
see that the Q-factor becomes large when most Rydberg 
atoms are paired up. Q ~ 1 already represents a siz- 
able deviation of the atom statistics from a Poissonian 
distribution. 



V. VARYING THE SPATIAL ATOMIC 
DISTRIBUTION 

The scenarios shown in the previous section are rela- 
tively weakly blockaded and assume a homogeneous dis- 
tribution of atoms. In this section we study how an in- 
homogeneous distribution and changes of the overall den- 
sity affect our results. The basic picture is shown in Fig. [4] 
and Fig. [51 In Fig. 2] we compare a variety of different 
densities. Prominent features of this data are compared 
more directly in Fig. [5j 

Is is known that the echo-signal, i.e. the fraction of Ry- 
dberg excitations remaining after the pulse, increases as 
the density is increased @. Using the w-expansion, we 
can supply the analytical expression Eq. for this in- 
crease. For the low densities where it is valid, it performs 
well as can be seen in Fig. [5] (a). Other effects of increas- 
ing density are a decrease of the maximum correlations 
and an increase of the Mandel-Q factor for intermediate 
times (r ~ 3/4). Despite the decrease of the correlation 
peak height, we find that the paired fraction R, defined 
by Eq. (fTTJ)) almost reaches unity after the pulse, regard- 
less of the density. This is shown in Fig. 2] (e). We also 
show the development of the dip in the correlation peak 
in Fig. [4] (f). We only see this feature for the lowest 
densities considered. 

The reduction of correlation strength with increasing 
density can be understood from the nature of the dis- 
tribution of excitations near the end of the pulse, using 
only classical arguments. Let us assume a homogenous 
distribution of pairs of excited atoms within a volume 
V, with a fixed distance ro between the partners of each 
pair. The distance rg then corresponds to the location of 
our pairing peak, ignoring its finite width. The orienta- 
tion of pairs in space shall be isotropic. Overall we thus 
have a distribution of positions for each pair: 



/(R,r) 



1 



4ttV 



%1-ro). 



(18) 



Here, R denotes the centre-of-mass position of a pair and 
r its relative co-ordinate. Assuming a total number of M 
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FIG. 4: (color online) Effects of inhomogeneities and changes of density on the pairing peak. We show densities poVbo = 5.4, 26.2, 
39 in (black, blue, green) respectively, and an inhomogeneous case p(x) = poexp [— 2r 2 /a 2 ] with poVbo = 26.2 and a/rbo = 1.6 
(red-dashed), (a) Rydberg fraction / e , for increasing homogeneous density the final remnant f e (r = 1) gets larger, (b) Paired 
excited fraction according to Eq. (|16[1 . (c) Mandel-Q parameter, see Eq. (|17[) . (d) Spatial maximum of pair correlations, (e) 
g^ 2 \r) at t = 0.75, for increasing homogeneous density, the peak height decreases, (f) Spatial atomic correlations near the 
end of an echo-pulse in a low density (poVi,o = 5.4) gas. Shown are snapshots at r = 0.96, r = 0.98, and r = 1.0 in order of 
increasing peak-height. 



pairs we have a pair density n p = M/V and an excited 
atom density n e — 2M /V. The positions of the pairs 
themselves are correlated due to the dipole blockade, we 
thus assume that the centers of the pairs must be sep- 
arated by more than a certain radius /■(,. The classical 
correlation function corresponding to our is 

9 (2 \r)= (19) 

f d n £, ( r i = x q; r 2 = x o + y> |y| = r ) 

J E{r 1 =x )S(r 2 =x + y, |y| = r) ' 

In the numerator we have the expectation value for the 
number of joint occupations of the locations n and r-i by 
excited atoms. In the denominator E(t\ = xo) denotes 
the expectation value for the number of excited atoms at 
location Xq, which is 2M/V. The integration dfl is over 
the solid angle of orientations of y with respect to Xo- 
We are now interested exclusively in the value of the cor- 
relation function at the pairing peak, located at r = ro- 
For simplicity we assume contributions to the numera- 
tor of Eq. (|19p stem only from cases where the atoms at 
both locations were constituents of the same pair. This 
should be justified when the pairs are sufficiently dilute 
that their mean distance is much larger than ro. We 
can then write E{r\ = x ;r 2 = x + y, |y| = r) = 
M/((r 1 +r 2 )/2,r 1 -r 2 ). 

With these simplifying assumptions, we find Eq. (119|) 
equals V/AM ~ l/n p for the simple pair distribution 
described above. Using n p ~ n e , n e = f e p and f e ~ p 



from the w-expansion, we find overall that the correlation 
strength scales like p~ 2 . This behavior is roughly con- 
firmed by the simulation results as shown in Fig.[5](b). It 
is consistent with a paired fraction R that shows almost 
no density dependence near the end of the pulse: For 
each excitation we calculate N = f 7 r ° +d r 2 g^ (r)n e (r)dr 
to obtain the number of excited neighbors in an inter- 
val [ro,r + d]. Since n e scales like p 2 , while scales 
like p~ 2 , we can understand how the final paired fraction 
can remain almost the same even though the density is 
varied. 

We verified that the relation g^ 2 \rQ, r) ~ p~ 2 remains 
qualitatively unchanged if we model similar echo-pulses 
in a fictitious system with a 1/r 4 long range interaction. 
This lends further support to the simple explanation in 
terms of the density of available pairs. 

Now consider a case with iVo = 50 atoms inho- 
mogeneously distributed with atomic density p(x) = 
po exp [— 2r 2 /a 2 } [case (iv) in tabled]. The peak-density 
po is chosen as for case (ii). Fig. @] includes the cor- 
relation function averaged over all atomic pairs in the 
cloud for this case. It is determined in the same manner 
as described in section IIIII One could expect that the 
inhomogeneity washes out the signal in the correlation 
function. Instead its visibility is even better than for the 
homogenous case with equal peak-density, owing to the 
presence of low density regions in the atomic cloud. For 
lower density the pairing effect is more prominent as we 
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FIG. 5: (color online) Dependence of echo signal and final cor- 
relation peak on atomic density, (a) Final Rydberg fraction at 
t — T. These datapoints are very well described by Eq. 
(b) (?) Peak height of g {2) at t = T. The (red-dashed) curved 
shows the functional form g^ 2 '(p) = [po / p] 2 g^ 2 \po) with ref- 
erence density po = 26.2. This is motivated by a simple model 
explained in the text. 



already argued. We also do not expect a strong spatial 
variation of the preferred distance throughout the cloud, 
since the blockade radius rb(t) depends only weakly on 
the atomic density (i.e. the saturated many-body block- 
ade radius scales like ft ~ p 2 J b Q @)- Hence the location 
of the correlation peak, situated near rb(r/2), is expected 
to vary only little throughout the sample. Consequently 
there are only small effects of spatial averaging on the fi- 
nal result in the inhomogenous case considered here. To 
fully exclude that density inhomogeneities in an exper- 
iment would suppress the signature reported here, one 
ideally should consider a sample whose width greatly ex- 
ceeds the characteristic range where the correlation peak 
is formed (a ro). This is however computationally 
intractable. 



VI. COMPARISON WITH ^-EXPANSION 



The simple expression Eq. (|T0[) for two-body correla- 
tions in the w-expansion does not depend on density and 
hence cannot capture its effect on correlation dynamics 
as shown in Fig. 2] However, for small densities it com- 
pares quite well with the substantially more involved ex- 
act Schrodinger evolution. This can be expected from 
Eq. (TT2"j) : For Nt = 40, the series is well described by its 
first two terms until ilt = 0.25. 

We can see in Fig. H (a) and Fig. [5] (a) that the NLO 
result for the Rydberg fraction [Eq. ©] gives good quan- 
titative results for the cases considered here. For corre- 
lations, we only have the LO expression Eq. (Ti"Tj|) . We 
see that this approximation describes correlations in the 
low density case poVbo = 5.4 quite well, while quantita- 
tive deviations appear, once the density becomes as high 
as in the case with poHo = 26.2. Near the very end of 
the pulse also qualitative differences in the shape of the 
correlation function are present, since Eq. (fT0| cannot 
describe the dip seen in Fig. 2] (f). 

We note that in particular for the cases that show 



the most dramatic correlation dynamics through an echo 
pulse, those with a low density, the w-expansion provides 
useful results. 



VII. CONCLUSIONS 



We have shown that an echo-type excitation sequence 
as employed in the experiments [§, 9] can be used to cre- 
ate Rydberg gases with non-standard nearest neighbor- 
distribution functions. After the pulse the vast majority 
of excited atoms possesses a neighbor in a fairly nar- 
row interval around some distance ro. Variations of the 
interaction-strength and pulse- length can control r . The 
strength of the correlation signal is proportional to p -2 , 
where p is the atomic density. This can be understood 
in terms of the quantum state after the echo-pulse, inde- 
pendent of the precise form of the interaction. 




FIG. 6: (color online) Comparison of the tj-expansion (red- 
dashed) with exact solutions of the SE for poVto = 5.4 (solid 
black) and poVto = 26.2 (solid blue), (a) Excited state frac- 
tions, (b) Pair correlation function g^ 2 \r) at r = 0.49. (d) 
Pair correlation function g 2 (r) at r = 0.6 and r = 0.75. The 
later time has a higher pairing peak. For r = 0.6 we used the 
same color scheme as in the other panels, for r = 0.75 the oj- 
expansion result is (magenta-dashed) and the exact solution 
is (solid green), (c) Spatial maximum of pair correlations. 



The described pairing effect in the density-density cor- 
relation function is most pronounced for low densities and 
not too strongly blockaded gases. However the overall 
fraction of Rydberg atoms that have a neighbor near the 
distance ro approaches unity regardless of density. Our 
results were obtained by direct simulation of ensembles 
of about a hundred atoms. Further we used these simu- 
lations to estimate the range of validity of the first order 
approximation of the correlation dynamics obtained us- 
ing the w-expansion. We find that for weakly blockaded 
gases it provides useful estimates. 
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